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1. Introduction 

The attempts to describe gravitational phenomena in terms of spatial gravitational 
fields are known from the times of Maxwell [T]. After establishing the General Theory 
of Relativity it became clear that the vector theory could not serve as an adequate 
description of gravitational interactions since the linear equations can be viewed at best 
as some approximation of an adequate theory, presented by the tensor equations of 
General Relativity. 

In this paper we suggest a modified version of a Maxwell-like system of equations for 
the gravitational field, calculate the defiection of a photon passing the central body and 
the precession of the planetary motion, and compare the results with those of General 
Relativity 0, El- 

We postulate fields g and B satisfying the system of equations 

dB 

div g = ArcGp, curl g = 

rfz.B=0, cnr/B=l|-^j (1) 

In (P), G is Newton's constant, c is the speed of light, p is the mass density and j is the 
mass density fiow. 

We also postulate that the gravitational force acting on a moving particle with 



velocity v and mass m (m = mg 7, 7 ^ = y 1 — v'^/c^) is described by: 

^ dp dim-v) n (va)v 

dt dt ' c2 

= mg + mvx(B + SB) (2) 
where a = dv/ dt, and 6B is defined by 

-(B.»)4f-l^J (3) 

Equations similar to (0) and to (without the 5B term) appear in the numerous 
attempts to model gravitational force after the electromagnetic Lorentz force (see e.g. 
Ill - El)- The system of linear equations (H)) for the effective spatial gravitational fields 
are not accounting for the self-interaction, which appears in the framework of General 
Relativity in a natural way, and which is essential for deriving the correct values for the 
observables such as the photon defiection and the perihelion precession 0-0. 

Introduction of the 5B-term via (0 effectively leads to an extra contribution to 
the B field which, in turn, results in an additional contribution mv x 6B to the force 
(everywhere below for notational simplicity we write B in place of B + 6B). 

Based on equations (0) - ©, we calculate to the lowest order in f^/c^ the defiection 
of a photon passing a central body and the precession of the planetary motion. It turns 
out that equations (jT)) - © yield results which coincide with those following from the 
General relativity |2], [2, suggesting that at least some observables can be calculated 
using fields g and B. 
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In section 2 we derive expressions for the angle of the deflection of a photon as 
it passes a central body (throughout, we assume that the central body is the sun and 
use a coordinate system in which the sun is at rest), and for the precession of the 
perihelion of a planet orbiting the sun. In section 3 we discuss our results as well as the 
limitations and approximations of the approach. In the Appendix we demonstrate how 
the equations of planetary motion which follow from our formulation and which lead to 
the correct expression for the precession, also follow from General Relativity provided 
that higher orders of v'^/c^ are neglected and characteristic distances are bigger than 
Schwarzschild radius. 



2. Deflection of the photon and the precession of the perihelion of a planet 

2.1. Photon Deflection 

We consider a photon of energy E moving along the y axis and passing the sun at a 
perpendicular distance d. At time instant t the angle from the x axis to the photon is 

e{ty. 

eit) = arctan ^, = ^ cos^^(t) (4) 

d dt d 

The x-component of the g field is 

= ~ COS a* = — COS (5) 

where M is the mass of the sun and R = d/ cos 6 is the distance from the photon to 
the sun. 

From the equations we obtain for the derivative of g 

The value of dg/dt is not zero because we evaluate the change in g along the particle's 
path; in this case the non-trivial time dependence is realized by 6' = 0{t). 

From the symmetry considerations it follows that the y component of the 
gravimagnetic field is vanishing: = 0, so integrating in (jH)) and using cos^ = 

d/{d'^ + y^)^/^ we obtain for B^: 

B. = / -^dy = - cos^'e (7) 



dy cd"^ 

Now, having g^^ and B^ at hand, we can evaluate the expression for F^. 

EGM . E GM . EGM „ 

F^ = mg^ + mVyB^ = —— cos^ 9 ^ c — -r- cos^ 9 = -2 — —— cos'' 9 (8) 

d'^ cd'^ d^ 

where we have used m = E/c^ for the mass of the photon, and Vy = c for its velocity. 

The coefficient 2 in (jSJ is due to the contribution from the "B" - term; as is evident 

from (jS)) the contributions from the "g" and the "B" terms are equal. In Newtonian 

physics Fa; would have been 1/2 of the value we obtained in (jH)), and since the deflection 

angle is linear in F^. ( see below Q), this would result in the half of the deflection angle. 
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We now calculate the change of the x component of the photon's momentum along 
its trajectory, 6px- 

roc AGME 

= / dtF, = — — (9) 

Therefore, for the value of the deflection angle 6 defined from the expression 
ta.n{6 (p) = \Spx/p\ "we obtain tan 6 if = AG M/c'^ d, i.e. for small 6 ip 

f4GM\ 4GM 
5^ = arctan(^-^j ^ (10) 

The same result follows from the General Relativity |2] , |S] • 
2.2. Precession of the Perihelion of Planetary Motion 

Let us consider the motion of a planet with mass m around the sun. At the moment t 
we have for the planet r{t) = {x, y, z) = (i?(^(t)) cos^(t), Ri9{t)) sin^(t), 0). 

For simplicity, we will present the calculations for the x - component of the force 
and the acceleration at the point 9 = 0. The results do not depend on this choice. 

From (0) and B^. = we obtain 

ox dt c^ 

where M is the mass of the central body, u = 1/R, and 6 = dO/dt. Using x = cos6/u 
and integrating (fTTj) over u leads to: 

B^ = ~^ue (12) 

Combining (0), JHl), and Vy = cosO/u - u' 6 sin6'/u^ = O/u (for 6* = 0; m' = du/dO) 
we obtain for the x - component: 

c^ 

Expanding in v"^ /(? and retaining the first non-trivial terms, from F^. = dp^/dt and 
(HSl) it follows that 

- GMmu^ - ^^me^ = (l + ^]ma^ (14) 
c^ \ c^ I 



where is the x component of the planet's acceleration. In Newtonian physics and 
Special relativity, the term mO'^ G Mj (? would not appear. 

Substituting = -u'e/u^ and = -u"e^/u^ - n'O/u^ + 2{u'ef/u^ - O^/u 
after straightforward algebraic manipulations we get 

~GM ^ + = i ^ + 2^^^ 15 

We can write expression (jlSp in a more compact form by using L = L/m = d/u^ where 
L is the angular momentum. We have 

GMu^U GMiu'Lf ,,-2 I ,dU 
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We will seek the solution of the equation of motion of the form u = 

uq + A cos{ke + B). 

From equations it follows that the y - component of the force (for 6' = 0) is 
non-vanishing, and is given by 

GMmu' ^2 

F^ = mB,v^ = 17 

We equate F^ to dpy/dt and neglect the higher order terms in 1/c^ to obtain for 



a. 



y 



dvy/dt: 



ay = 2 18 

Note, that in the absence of the B field we would have Fy = 0, and as a consequence, 
ay would have been half as large. 

Combining (|T8|l . ay = udL/dt, and ansatz u = uq + A cos{k9 + B) we obtain 

for dP/de 

dU GMu' GMPA, ^ , 

4 - = 4 ksm{k9 + B) (19) 



de * c2 c2 
In approximation GMA/c^ ^ 1 the solution of (fT^ is 

= Lo' (^1 - cos(A;^ + B)^ (20) 

Now we substitute u = uq + A cos{k9 + B) and (j20|) into equation (fT6|) . and equate 
the same powers of cos{k6 + B). The solutions for uq and k"^ are as follows: 

Uq ~ — ^, (21) 
-^0 



and 



e^i-Q^ (22) 

^0 



The corrections to ()21|) - ()22|) are of higher order in 1/c^, and are numerically negligible. 

Expression (|22j) for k"^ coincides with that from General Relativity [2], [S], and 
of course, yields the correct result for the precession of the planetary perihelion. In 
particular, expression ()22p gives the correct value for the observable precession of the 
perihelion of Mercury. 



3. Discussion 



We have demonstrated that the postulated approach, equations (HJ - ©, leads to the 
correct expressions for the two key classic tests of General Relativity - the precession of 
planetary motion and the deflection of a photon. 

The B field leads to a doubling of the force on the photon, and therefore to 
a doubling of the deflection obtained in the framework of Newtonian gravity. This 
doubling leads to the correct result for the deflection. 
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In the case of planetary motion, the B field leads to an additional term in F^. and 
to a doubling of ay (see Equations (fTSjl and (fTHj) ). These values of F^. and ay, used in 
the equation of motion (jlfjj) . lead to the correct expression for the perihelion precession. 
If the B field were ignored in calculating F^. and ay, we would have obtained for k'^ the 
expression 1 — 2G^M^/c^Lo^ instead of the correct value /c^ = 1 — 6 M^/c^ Lq^- 
The observation that equations (HJ - © yield the correct results for these two different 
phenomena seems more than mere coincidence. 

Let us address an important question about our approach: what is the justification 
for the key postulates (HJ - ©• In addition to the observation that they lead to the 
correct expressions for observables, we note that equations fllH|l and (fTTj) which are 
immediate results of our approach (Q) - Q follow from General Relativity, assuming 
spherical symmetry and time independence of the metric. This derivation demonstrates 
consistency of our approach with General Relativity in case of weak fields and in order 
jc?. We sketch this derivation in the Appendix. 

As an attempt to interpret our approach, we note that switching from a partial 
derivative to a total derivative in (jS)) which introduces an extra B term and which leads 
to the correct expressions, can be viewed as calculating the changes in the fields as 
experienced by the particle, rather than the changes occurring at a point. It is plausible 
and reasonable to at least speculate whether the force acting on a particle might be 
related to the fields and changes in those fields as experienced by the particle. 

The approximations we used refiect a requirement that the gravitation field is weak. 
For example, from our expression for the photon defiection, the requirement that the 
corrections be negligible implies that 



When this condition is satisfied the corrections to our results ()10|) and (j^^ are negligible 
and we are able to reproduce the results of General Relativity. In terms of the length 
scales the above requirement reads that 



where Rs = 2G M/c^ is the Schwarzschild radius for a source with mass M. When ()24j) 
is not satisfied, our results differ from those of General Relativity; e.g. the expansion of 
arctan 4 G Mj c? d in (fTUI) generates the correct functional form for the next correction 
photon deflection angle, namely {AG M /c^ dY , but with the wrong coefficient —1/3. 
The same is true for the precession of the perihelion. This higher-order discrepancy 
is not surprising since our formulation is not proposed as a replacement for, or an 
improvement of General Relativity. Rather it should be viewed as an attempt to describe 
the gravitational phenomena using Maxwell-like system of equations. An advantage of 
the suggested formulation is that in a leading approximation it reproduces the correct 
expressions for observables using simple equations Q - (0); without invoking tensor 
calculus. 




(23) 



(24) 
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Though the results (fTUI) and (j221) were obtained from General Relativity years ago, 
it is quite unexpected that they can be obtained using an approach other than one based 
on Einstein's equations. This alternative "heuristic" derivation, and not the expressions 
()10|) and ()22|) themselves, is what we consider as the main result of our paper. From 
our viewpoint, the possibility of an alternative derivation of expressions for observables 
such as photon deflection and perihelion advance may provide an impetus for further 
understanding and insights about the theory of spatial gravitational fields and General 
Relativity. 

To conclude, from the equations ((T)) - Q it is possible to obtain the correct results 
(in leading approximation in Rs/R) for the photon deflection angle and the precession 
of the periastron. 
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Appendix. Derivation of Equations (|13|) and (|17p from General Relativity 

Let us sketch how eqs (fTT^ and (fT7j) can be derived from the equations of General 
Relativity. These equations also follow from our formulation, and lead to the correct 
expression for the precession of the perihelion. 

We start from the isotropic metric expressed in geometrized units as 

.oo--(l-^ + ^j, ..-^.(l-^j (25) 

Metric fl25|l is obtained by applying coordinate transformation i? — » r (1 — 3M/4r) to 
a Schwarzschild solution g^y{r). In the order of approximation used, terms of higher 
order in gjj are negligible. 

The motion on geodesies is described by 

+ KpPy^ = (26) 
where is the a component of the momentum, p*^ = dp''/dx°', and F^^ are the 
Christoffel symbols 013. 

For the x - component we have 

"^0^ + r>V = 

where r is affine parameter of the geodesic. The non vanishing Christoffel symbols are 

F^ -d-^]^ F^ --(l-^JiV^ 
°° " 1^ 2r) i?2' - R J 

i.e. 

Neglecting higher order terms of M/R we have (p^)'^ ~ niQ (1 + 2M/R + p^/mg) which 
when substituted into ^T7\ leads to the equation of motion in the xy plane: 

Using p = mv and neglecting higher order terms of v'^/c^ leads to 

and, using dr = y/—goQ dt ^ 1 — M/R, 

dp- moM A ^ M\ moM ^ .^n^ 

Recognizing that 7^1 + M/2R, R = u^^, m = '-fniQ, and v"^ = u^^9 (at = 0), we see 
that equation (jHUjl is the same as our equation (fT!?jl . 

Equation (fT7|) is obtained in the same manner, considering the y component of 
equation (j2l 



